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ABSTRACT 


A discussion is given of particle lags in mean flows, acoustic 
oscillations at single frequencies and in turbulent flows. Some 
simplified cases lead to exact solutions. For turbulent flows lineariza- 
tion of the equation of motion after assuming the fluid and particle 
streamlines coincide also leads to a solution. The results show that 
particle lags are a function of particle size and frequency of oscillation. 
Additional studies are necessary to evaluate the effect of turbulence 
when a major portion of the energy is concentrated in small eddies. 



INTRODUCTION 


A promising method of measuring turbulent flow characteristics is 
the Laser-Doppler technique. Local mean and fluctuating velocities of 
small particles suspended in a fluid stream are determined. The seeding 
of small particles is necessary to provide scattered light for the Laser- 
Doppler instrument. The question then arises: "How well do the particles 
follow the flow?" This report considers the motion of a single particle 
or group of particles in a fluid stream with emphasis on solid particles 
in a gas. Flow in pipes and in an isothermal round free jet are used as 
examples . 

Motion of a single particle is a simplified case of dilute gas-solid 
suspension flow. Large particles will only follow the slow large scale 
turbulent motions of fluid, while particles small compared to the smallest 
scale of turbulence will respond to all turbulence components of the fluid. 
Since the particle contributes to energy dissipation due to lags between 
particle and fluid, large numbers of particles will modify the energy 
spectrum of the fluid. 

Many studies on particle fluid interactions have been made because 
this problem is of importance in several fields. The motions of small 
particles have been analyzed by comparisons of the amplitudes of oscilla- 
tions of particle and fluid in an accoustic field, and by comparison of 
the mean velocities in turbulent flow. Special problems exist when solid 
boundaries are present or particle-particle interactions are important. 

At speeds below the velocity of sound the difference velocity between 
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fluid and solid is near the continuum regime and viscous effects are 
important. Corrections for compressibility and rarefaction are necessary 
for higher speeds when slip or transition flows are encountered. 

The results of calculations on particle fluid interactions for small 
particles generally show that individual particles below one micron in 
diameter follow the flow. Turbulence measurements on particles are then 
a good representation of the fluid behavior as long as the turbulent 
fluctuations are not extremely rapid. Further work is necessary to 
consider the energy dissipation range of some turbulent flows or high speed 
f 1 ows . 

In the following pages the velocity of the fluid and particle are 
compared for several conditions of physical interest. First the mean 
velocity and then the velocity fluctuations are calculated for a single 
particle. A survey of the literature on experimental work is also pre- 
sented. Finally some consideration of multiparticle systems is given. 
Although the equations are the same, the conclusions for flow in a liquid 
are much different from those given in this report for a gase. 

Extensive introductions to the field are available in the books by 
Fuchs [1] and Soo [2]. These references have been used extensively in 
this report. 



I. ONE DIMENSIONAL MEAN VELOCITY OF SINGLE PARTICLES 


We begin with the simplified case of one dimensional flow when the 
particle does not distrub the fluid motion. 

The equation of motion for a single spherical particle in creeping 
flow is given in many texts as: 
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dV 

dt 


= 6TT U r (U-V) - ^ & 


p dz 


nr 
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d V _ dU 
dt ” dt 


2 

6r /irup 


d(V-U) 

dx 


dx 

/t^ 


(1.1) 

+ F_ 


This equation was originally obtained by Basset [3] and is discussed by 

Tchen [4] and Hinze [5]. The notation is as follows: 

4 3 

m = mass of particle = j irr Pp 

4 3 

m' = mass of fluid = j Trr p 
V = particle velocity 
U = fluid velocity 
y = fluid viscosity 
p = fluid density 
Pp = particle density 
r = particle radius 
t = time 

• F g = external potential force 

Equation 1.1 then equates the force to accelerate a particle to the follow 
ing in order of their appearance in the equation: 

1. Viscous drag on the particle for slow relative motion. 
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2. Pressure gradient in the fluid surrounding the particle. 

3. Force to accelerate the apparent mass of the particle relative 
to the fluid. 

4. The effect of acceleration on the viscous drag. (Basset force) 

5. The external potential force such as gravity. 

To account for instantaneous motion in a turbulent field, this 
equation must be written in vector form in the proper frame of reference. 
However, the equation as it stands can be used to determine the behavior 
of the mean flow. In a later section the particle velocity in a fluctuat- 
ing medium will be considered. In general the pressure term can be 
obtained from the equation of motion of the fluid assuming the particles 
do not have an effect on the fluid behavior. The equation then depends 
upon the viscous forces and inertial forces in the fluid and the results 

will be functions of velocity gradient (dv/dz) and the gradient of velocity 
2 2 

gradient (d v/dz ). In this discussion these terms are assumed small or 
that any effect can be included in an experimental drag coefficient. 

That is the viscous drag term in Equation 1.1 is replaced by mF D (U-V) 

where 

f = ! c d^< u - v > 

and C D is the drag coefficient. At low particle Reynolds number, a 

the first term in Equation 1.1 is obtained. 

If the particle is in a fluid moving with constant velocity, the 
Basset force is neglected and the pressure gradient and external potential 
force are constant, the solution to Equation 1.1 is: 
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U-V 


2r 2 
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_ 9yt 

r 2 (2p +p) 
e H 




' ( 1 . 2 ) 


when at t = 0, the particle and fluid have the same velocities. The largest 

2r 2 


difference then occurs at t = » and the factor 


9y 


m 

dz 


need only be con- 


sidered in the absence of a potential force. For flow in a pipe the pressure 
drop is a function of the pipe diameter, D, and average velocity, <v>. In 
terms of the Fanning friction factor, f, 

2f 


-£*£ <*>• 


The friction factor is given by: 

f = 0.0791 N Re % 

at Reynolds numbers between 10 4 and 10 5 but for purposes of comparison it 
can be assumed that the result holds for higher fluid Reynolds numbers. 
Then: 


(U - V) 


00 


< V > 


.0088 




(1.3) 


where R is the pipe radius. To obtain less than 1% relative error at 
N Re = 10° in a 4 inch diameter pipe a particle less than 300 microns in 
radius is required. If the pipe is vertical, the potential force due to 
gravity must be included. 



F e . = (P p " p) 9 


6. 


Then: 


(U-V)c 

<v> 


, , o r y(p D ~ p) 

- - 0088 <*> V 7 ' * |r— ^ — 


-r N 2 




2r 2 gp p 

9y<v> 




o _4 

When Pp = 1.0 g/cm and y = 2 x 10 g/cm sec., 

10 6 r 2 

IT ^ 

The gravity force is larger than the pressure force and results in a lower 
particle size limit for a 1% error. This limit is approximately equal to 
<v> 2 microns when the velocity is given in cm. /sec. For example under the 
same conditions as above a 95y diameter particle is the 1% limit while a 
300 micron particle would give an error over 10% when the fluid stream is 

O 

air and the particle density is 1.0 g/cm' 2 . 

The analysis appears valid for pipe flow mean velocities in isothermal 
turbulent flow when compared with experimental data. Recently Reddy, Wijk 
and Pei measured particle and fluid velocities in a 10 cm. diameter vertical 
pipe [6]. Using 200 micron particles the error at the centerline was 16%. 
Assuming they used alumina particles with a density of 2.42 g/cm 3 the cal- 
culated error at t = <® is 25%. The agreement is good and could be improved 
by using a correction to the viscous drag. 

If the terminal velocity is not reached, the particles may show varia- 
tions in velocity due to different times of approach to terminal conditions. 


o 
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It is important to keep the velocity error less than 1% if this is to be 
avoided. 

Hughes and Gilliland [7] showed the results of Basset's solution to 
Equation 1.1 for motion of a particle in a fluid at rest. The terminal 
velocity is not changed by the addition of the integral term (Basset force) 
when the pressure gradient and external forces are constant and the fluid 
velocity is also constant. This can be easily shown by the Fourier trans- 
form method which is explained in Chapter II of this report. The tran- 
sient curves, however, are changed. The plot given by Hughes and Gilliland 

3 

shows little difference in the transient curves when p p /p ^ 10 but for 
Pp/p < 10 the corrections are significant. 

The next problem to be considered is to add acceleration of the 
fluid as would be encountered in the free jet. Again the Basset term will 
be neglected. Also the external force will be dropped and the assumption 
made that: 

dp _ dU 
“ ciz p dt 


from the Navier-Stokes equation of fluid motion. Then: 


| r 3 (2 P ♦ p) 


dV _ dU 
dt ~ dt 


+ 6yr (V-U) = - j r 3 (p - p) ^ (1.4) 


If the velocity varies linearly with downstream distance, Z, 


U = U 0 + 6z, 


(1.5) 


where 6 is a constant. 
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Since 

U = dz/dt 

< 

3T= «« ' 0-6) 

And 

U = U 0 e 6t (1.7) 


Solving Equation 1.4 then gives 


U-V 

U 


26( p p - p ) 

(2p p + p)[e + 5) 


[1 - e 


-(6+e)t-j 


Where 


e = 


9u 

r 2 (2p p + p) 


( 1 . 8 ) 


At the centerline of a round free jet the mean velocity decreases almost 
linearly for the first 21 diameters after the velocity decrease begins as 
shown by Corrsin [8]. In this region the parameter <5 becomes approximately 

Uo 

6 = - 0.05 -g— , 


where U 0 is the velocity at the jet exit and d the jet diameter. Errors 
are less than 1% for particles less than 10 microns when the particle 

3 

density is one gram per cm for <5 = 10. Similar calculations have been 
made by Gilbert, Davis and Altman [9] for the transient case. In super- 
sonic flow 6 > 500 and particles less than one micron radius are necessary 
to have less than 1% error considering an infinite time. 
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Inclusion of the Basset force will effect the transient behavior 
but not significantly the terminal velocity difference at infinite time. 
Measurements of particle mean velocities will, therefore, closely approx- 
imate the mean stream velocity when the particle size is on the order of 
one micron. 

A general solution to equation 1.4 including the gravity term was 
given by Tchen but the solution is difficult to apply. Most measurements 
are taken under conditions where the time from the flow source is much 
shorter than that to reach the terminal velocity. Therefore, errors in 
mean velocities are less than the figures used in the illustrations above. 

Deviations from the creeping flow assumption occur when the particle 
Reynolds number is greater than 0.1. Empirical equations can be used to 
modify a drag coefficient in Equation 1.1 for these higher Reynolds numbers. 
The results will lead to a reduction in velocity lag and are discussed in 
the section on turbulent flow. Numberous theoretical studies are avail- 
able in the literature on the flow of a compressible gas containing solid 
particles. A review up to 1966 has been given by James, Babcock and 
Seifert [10]. In this report the effects of compressibility and rarefac- 
tion will only be considered in the section on turbulent fluctuations. 

Only a limited number of experimental studies have been made in 
which the particle and fluid velocities are measured. These are listed 
below with a short description of the results. 

1. Recent Studies of Mean Velocities in Pipes 

1969 - Reddy, Van Wijk and Pei [6] - Measured air and 



particle velocities in 10 cm. tube at Air Reynolds 
number of 55,000. At centerline 200y particle velocity 
was 8.13 m/sec compared to 9.8 m/sec for the air. 

1968 - McCarthy and Olsen [11] - Summarizes work on 
drag coefficients and measures particle and fluid 
velocity in a pipe for 64 micron glass beads and 3 
micron CaC0 3 in accelerating flow. After the initial 
entry the velocity of the 3y particles was indistin- 
guishable from the gas while the 64y particles had a 
velocity 10% less than the air velocity. Using the 
formula 

dv 3C D p(U-V) 2 
dz 8p p rV 

They found good agreement with the experiment when Cp 
was the standard drag coefficient. Their results 
support the work of Torobin and Gauven [12] who found 
that wake interaction and free stream turbulence can 
be neglected in dilute particle suspensions. 

Recent Studies of Mean Velocities in Nozzles or Jets. 

1969 - Morse, Tull is, Seifert and Babcock [13] - Used a 
laser Doppler method for nozzle flows. In their A1 2 0 3 
system errors of 6.8% for a one micron particle to 
41.9% for a 7 micron particle were predicted. From the 


particle size distribution they matched the signal 
amplitude vs velocity. However, signals at lower 
velocities were unexplained by a collisionless 
model. Computations based on the effect of colli- 
sions showed that smaller particles would undergo 
large velocity changes in collisions with larger 
ones. They concluded that particle size distri- 
bution must be known in order to calculate the 
velocity distribution and particle collisions must 
be included in the model for calculating velocity 
lags. 

1968 - James, Babcock and Seifert [14] - Used a 
laser-Doppler technique to measure nozzle flows. 
Indicated that particle collisions are important 
at low velocities. Results agreed with calculated 
according to the authors but no quantative values 
are given. 

1965 - Fulmer and Wirtz [15] - Measured photograph- 
ically the particle and gas flow at a nozzle exit 
as a function of particle size. The gas velocity was 
2,900 feet per second. Figure 1.1 is a plot of the 
data for the aluminum powder velocity - gas velocity 
ratio vs the particle diameter. 
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Other references are listed by Soo [2] for the pipe flow case but no others 
could be found for the nozzle flow. Measurements on turbulent spectra and 
fluctuations will be discussed in a later section. 

Another factor often questioned is the effect of particle rotation. 
Afanasev and Nikolaevskii [16] have shown that the rotational oscillations 
of a solid particle decay exponentially with characteristic relaxation 
time 

= i_i£ 

T w 15 y ^p 

Small particles in air then will equilibrate very rapidly to the condition. 

3 

of the surrounding fluid. A one micron particle with density of 1 g/cm 
has a decay proportional to exp(- 3 x 10 5 t) when t is in seconds. Only 
particles over 100 microns in radius would maintain their rotation for a 
sufficient time to effect results unless very short times are investigated. 
In mean velocity measurements, rotation is not significant; but for tur- 
bulent fluctuation measurements these decay times may be significant in 
small eddies. 

Rotation of a solid particle in a fluid stream can be created in 
many ways including: 

1. Velocity gradients in the flow field. 

2. Collisions with solid boundaries. 

3. Collisions with other particles. 

A rotating particle adds a lift force to the equation of motion but 
does not effect the drag coefficient unless the relative Reynolds number 
is high. At the high Reynolds numbers rotation thickens the boundary layer 
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and changes the characteristics of flow separation. [17] 

Some studies of the effects of rotation have been made by Rubinow 
and Keller [18], Eichhorn and Small [19] and Saffman [20]. Although 
Saffman concludes the result of Rubinow and Keller is high, their lift 
force is to the first order 

F l = trr 3 p^ x H 

where is the vector angular velocity. For small spheres on the order 
of one micron in diameter in air, the angular velocity must be very large 
before this term is significant. 

The ratio of the magnitude of lift force to drag force is approxi- 
mately [2] 

| | _ 1 r^wp 

| f d| " 6 k ‘ 

For a one micron particle in air at a rotational frequency of 20 KHz this 
amounts to only 6%, but for particles two microns or larger this force is 
significant at high frequencies. Another discussion of the left problem 
has been recently given by Kondic [21], 



II. MOTION OF PARTICLES IN AN OSCILLATING FLUID 




Although turbulence is composed of random fluctuations, before an 
attempt is made to consider turbulent fluctuations, it is well to examine 
fluctuations of a single frequency. The one dimensional equation of 
motion of a particle in a vibrating medium consists of Equation 1.1 plus 
terms due to the oscillation. Additional terms, one proportional to the 
velocity which increases the drag and one proportional to the acceleration 
which increases the inertia part are required. The extra forces added to 
Equation 1.1 are: 


F osc = - 




(V-U) - f m'B 


dV _ dU 
d"t Vdt 


( 2 . 1 ) 


Where is the frequency in sec. -1 and 8 = -p /2]i/po> 0 . This formula is 
derived in Lamb [22], and Landau and Lifshitz [23] based on the earlier 
work of Stokes. The equation of motion is then: 


m = m ' <31 " f m' u) 0 8(l + 8)(V-U) 


- m' 


1+ i b 
2 4 3 


d V _ dU 
dt " dt 


- 6r /irpp 


d(V-U) dx 
dT /XX 


( 2 . 2 ) 


or in terms of A = U-V 



1 9 \ 

m + m' 

{ zU\ 


^ j m'u) 0 8(l + 3)A + 6r 2 /rryp 


I 1 


dA dx 
dT 


( dU 
(m-m ) ^ 


(2.3) 
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This equation is easily solved using a Fourier integral transform: 

r 


F(o>) = 


e“io)tf(t)dt 


‘(2.4) 


and the inverse transform 


f <‘> - k 


e 1a) ^F(co)dai 


(2.5) 


We assume U = U Q sin w t and t Q , the starting time is minus infinity. 

Then the integral term has a simple transform and the transform of the 

sin oj t and cos <a t are Dirac delta functions. The inverse transform 
o o 

needs no integration. F((o) is evaluated at two points to give f(t), or 
the amplitude of the response can be obtained directly from the square 
root of the product of the transform with its complex conjugate. 

Let us write Equation2.3 as: 
rt 


dA . f K , 
£ dt + ^ + T 


= U w„(m - m‘ ) cos w t 

dx rz — o o' o 

/t-T 


( 2 . 6 ) 


If A = 0 at the starting time (-«>), the transform, A(u) is 


^ j Vo (m " m ') 7r ^ (S ( 0J " w 0 ) + + »o>3 

iwe + £ + y(l + i)w^ /tt/2 


(2.7) 


Then the inverse after some manipulation becomes 


a) (m - iii') 
o 

u + ^ Q h ^m 2 


— x 

+ (we + 6a) Q ^ /tt/2) 2 
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( 2 . 8 ) 


£ + Soo^ /tt/2 cos co Q t + Ioo£ + 6co 


O^) 


sin w t 


Or: 


A(t) _ 


co (m - m') 
o ' 


(2.9) 


0 £(£ + 6oo Q ^ /tt/ 2) ^ + (coe + 6co 0 ^ A/2)^J^ 


sin (co n t + c|>) 


Where: 


^ + 6ui 2 A/ 2 
tan * = 0 


coe + 6 co q ^ A/2 


( 2 . 10 ) 


In terms of r, p, p p and y the amplitude is then: 


( 2 . 11 ) 
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u U f 
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which is a function of p p /p and 8 only. This particular form of the 
comparision is not found in the literature and it is convenient to solve 
Equation 2.8 for V(t) to compare amplitudes directly. The result is 


-jj = V Q sin (cot - ^) 


( 2 . 12 ) 


Where 


V o = 


/ 2(^+l+f8 + — 


1 + A. 


1 - 


4/2- 


8) 




(2.13) 


and 


tan \p = 


(2.14) 



\ + 
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To compare with previous results we consider simplifications of Equation 
2.2 and their solutions. 

If all the terms on the right hand side of Equation 2.2 are dis- 
9 2 

carded except ^ m'ooB (V-U), the resulting equation represents motion in 

a medium whose resistance is non-inertial . Then for U = U„ sin wt, 

o 

m cHF + T m * m‘ co$^U 0 sin wt (2.15) 


The amplitude of the solution is 



(2.16) 


Equation 2.16 was used by Rosenweig, Hottle and Williams [24] to justify 
their light scattering measurements of smoke to determine turbulence 
fluctuations. Fuchs [1] cites several experimental confirmations of this 
expression and also shows that the solution to Equation. 2. 2 without the 
integral term gives almost identical results. Since the experimental 
data scatter widely, there is some doubt of the validity of the analysis. 


The solution to eq. 2.2 without the integral term has been widely used 
since its derivation by Konig [25]. The result is 



1 + 3g + 9 g 1 2 + £ b 3 + i g 4 ** 
_g 2 + 3gB + | B 2 + | B 3 + | 6 4 _ 


(2.17) 


1 

where g = -j (2-jj + 1). Since in all cases the amplitudes are functions 

of r^ET and pp/p, plots can be made by varying these two parameters. Figures 

2.1 and 2.2 show such plots. 
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In the first, Figure 2.1, Equations 2.16 and 2.17 are plotted 

3 

for air at room temperature and for a particle density of l.Og/cm . 

Figure 2.2 shows equations 2. 12, and 2.13 for the same conditions. Equa- 
tion 2.16 would apply only for low relative Reynolds numbers between 
particle and fluid and in the absence of other forces. Fuchs lists additional 
forces such as hydrodynamic sound pressure due to changes in density of 
the gas which cause particles in a standing wave to drive to the nodes. 

Other effects are scattering and absorption of sound waves, hydrodynamic 
interaction between particles and electrostatic scattering of particles. 
Little is known about the importance of these forces as functions of 
concentration, mean velocity or frequency. 



III. MOTION OF A SINGLE PARTICLE IN A TURBULENT FLUID 


In turbulent flow the oscillations do not occur at one frequency 
or amplitude but vary in time and space. Although the acoustic wave gives 
apparent agreement with single frequency experiments, it would not account 
for the frequency spectrum found in turbulence. If the measuring instrument 
is able to determine precisely the instantaneous velocities, the only errors 
in using the results for the fluid velocity will be the difference between 
the fluid and particle. Most instruments measure the velocities over a 
finite volume and a finite time span. Therefore a distribution of veloc- 
ities is obtained. The problem of investigating the relationship between 
the measured distribution and the true distribution involves the simul- 
taneous consideration of the instrument, the particle lags and the fluctu- 
ations due to turbulence. Only the interaction between the particle lag 
and the turbulent fluctuations are considered in this report. A possible 
way to allow for turbulence is to assume that Equation 2.17 also applies 
to the ratio of amplitudes of the frequency spectra curves of particle 
and fluid turbulence. 



Then 
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where F(oj) is the amplitude of the normalized frequency spectra density 
curves for the fluid and V 2 ", and U 2 " are time mean square averages. A 

i 

similar expression can be obtained for the relative velocity = V-U. 
We will compare Equation 3.1 with another analysis for examples of F (oj) 
obtained from experiments later. 

As in the ease of the acoustic field the equation of motion of a 
particle in a turbulent field is obtained from the momentum balance. 


where 




Forces , 


dV aj 
dt at 


(v.v)v 


in vector notation. For simplicity the vector symbol will be dropped. 
For a spherical particle when continuum equations are valid the forces 
are the following: 

1) Gravity = mg 

2) Buoyancy = m'g 

3) The variation of hydrostatic pressure in the fluid = m' — 

dt 1 

where the substantial derivative is based on the fluid 

velocity (indicated by — = — + U.v) 

dt 1 at 


4) The force necessary to accelerate the "apparent" mass of the 


1 


particle relative to ambient flow = 


m 


dU 

dt 


dV 

dt 


5) The deviation of the flow pattern from steady state (Basset 
force) = 

dU _ dV 
dr dx 

Vt -T 




24 . 

6) The drag force = mF D (U-V) 

7) The Lift due to spin. 

8) The pressure = ^-rrr^vp 

9) Intermolecular forces. 

10) Forces induced in the neighborhood of an interface. 

11) Forces due to temperature gradients in fluid or particle. 

12) Forces due to molecular diffusion. 

13) The force due to oscillating pressure 

14) Other body forces 

Only 3, 4, 5, 6, and 8 are usually retained in the equation of motion. 

These terms form an equation which has been solved only when additional 
assumptions are valid. Here we take a modification of the approach given 
by Soo [2] based on solutions by Tchen and Chao. The assumptions are 

a) The particle is small and spherical so that a simple 
expression can be used for the drag coefficient based 
on the relative motion. The drag coefficient appears in 

6) above when F Q = | |U-V| . 

P 

1 » ' n ar.Torobin and Gauvin [26] have reviewed the effect of tur- 
bulence on Cp. Literature results show both increased 
and decreased drag coefficients. An empirical correction 

• which fits the standard drag curve over a large range of 
particle Reynolds numbers is given by Torobin and Gauvin 
and is convenient for computer calculations. 

c ° = 4r (1 + 0,15 NRe °* 687) (3,2) 
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where Np 2ruak . 

Re y 

At high speeds approaching sonic velocity the drag co- 
efficient can be modified to include effects of gas 
compressibility and rarefaction [27]. If the effects 
of velocity gradient and derivatives of velocity grad- 
ient can be expressed as a constant multiplied by the 
relative velocity, the drag coefficient can be further 
modified. 

Cp (modified) 

Let K = C D (Stokes) 5 


Then only the factor K will appear in the equation of 
motion. K is a function of the variables mentioned 
above. 

b) The particle is small compared to the smallest wavelength 
of turbulence. 

c) The flow field is not perturbed by the presence of the solid 
particle. Then the components of the pressure term can be 
taken from the Navier - Stokes equation for the single phase 
fluid. 


where 



mV 2 U, 


DU _ aU 
Dt ” at 


U-vU. 


When a cloud of particles is present this equation includes 
the effects of flow around the solid. A modification of this 
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can be made if we can include the effect of the nonlinear 
terms in the drag, coefficient. 

d) The path of the solid particle coincides with the fluid 
streamline. That is = ^jp- . 

This assumption leads to the conclusion that the particle 
diffusivity is the same as the Lagrangion eddy diffusivity 
of turbulence. Soo discusses the experimental and theoret- 
ical reasons why this is not true for a finite size particle. 
The particle diffusivity is a function of particle response 
time and the Lagrangion and Eulerian scales of turbulence. 
However, without the assumption that the particle and fluid 
paths coincide, no rigorous solution to the equation of 
motion is available [28]. 

e) The turbulence is homogeneous and steady. Then a simple 
velocity spectrum can be used. Also the mean flow can be 
omitted from the equations by having the coordinates follow 
the mean motion. This assumption has been discussed in 
Section I. If we use a coordinate system which moves with 
the particle, the total time derivatives of the fluid velocity 
must be replaced by the derivative moving with the coordinate 
system: 

ft- • $ ♦ (V-»)U <3 ' 4) 

f) The turbulence is not effected- by solid boundaries and the 
starting time, t Q , in the Basset force term can be set at 
minus infinity. 
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With these assumptions the equation of motion becomes 
dV = 67 ryrK(U-V) + m' [|j + (U-v) U - W 2 U] 


m 


dt 


i Hit h\/ o r* d(U-V) dx 

+ & - & ♦ 6 ^/ 


(3.5) 


Or 


dV 

dt 




2 Pp+ P r 2 


}2p p+ pj 2p P +pdt 


t £ (U-V) 


2p p +P| r 


r~ C d? 

ftZ 


dx 


(3.6) 


Now .let a = > and 6 ~ 2 3 + » then: 

pr P p p 


t d 


4|r - 1^- Vp + a6K(U-V) + 


MU 6 /E [ Sr 


(U-v) 


aX V 7 T J 

—00 


dr 


/t- 


(3.6a) 


To write the components of equation 3.6 in a cartesian coordinate system 
we have 


dV. 3V* jii 

— 1 = — + (V*V) V. = (^) 

dt at i dt i 


(3.7) 


and we note that 


f£ + (u-v)u 


fjr + (V-V)u + [(U-V)’V]U = ^r + [(U-V) *V]U. 


(3.8) 
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So that 


[ an 1 dU -i dU i dU i 

3U + (U*v)U| = — + (U: - V.) — '- + (U. - V.) — ^ + (U. - V.) — 

3t Ji dt 11 3X- J J ax. k 3X 

J 


9U, 


3U, 


Then equation 3.6a becomes 


dV. dU. 9 3U, 9 3U. 

— - = 6 — - + a 6 K (U.-V.) + % 6 (U ,-V-)— U f 6 (U.-V.) — 
dt dt »y J J J av 


3X i 


3X. 

J 


2 3U i 2 2 3„ ^ d7 ( U i- V i) dT 

+ ! s < u k- v k> ^ % + s ^ 2 J ^ - < 3 - 9 > 


This derivation is also given by Hinze [5], and it represents the 
particle motion from a Lagranian point of view. That is the particle 
follows the detailed motion of the fluid. Friedlander [29] suggested a 
Eulerian formulation when the particle is large. The true situation lies 
in between the Lograngian and Eulerian models. A short discussion of this 
point and a review of the literature pointing out the important parts of 
the preceeding derivation is given by Brodkey [30], 

Equation 3.9 contains three nonlinear terms and the term involving 
the second derivative with respect to the space coordinates. Conditions 
for the neglect of these terms have been given by Corrsin and Lumley [31] 
and Hinze [5]. A similar argument is given here. Let us assume that 
(U.-V.) and (U.-V. ) are porportional to (U.-V.). Then the nonlinear terms 

J J K K I 1 

can be grouped together 



I'fitUrV 


9U i’ 


311. 3U. 

— + a. — + a, — 

. 9X i ' 9X j 9X kJ 


-T 6 < u rV A 
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(3.10) 


where a i and a^ are constants. 

If a6K(U^-V^) is much larger than this term, we can neglect the non- 
linear term. The condition becomes if aj and a^ are equal to one and all 
the derivatives are equal 
3U. 

26(u.-V.) —1 « a6K(U.-V.) 

1 1 3X i 11 

or since K is near one, . 


3U. 

2 — « a 


(3.11) 


Since a is approximately 10 7 for a one micron particle in air, the velocity 
gradient must be extremely large to be significant. Of course for a lOOy 
particle or a somewhat smaller particle in water this term may be important 
It is possible then to calculate a particle velocity error in the same 
manner as shown below by taking the local value of A and using it to form 
a new correction K. 


K 11 = K + ~ 
3a 


(3.12) 


This approach has been taken by Tchen [4]. The second derivative term can 
be neglected if it is small compared to the others. Hinze compares it to 
the velocity gradient term which gives 
V 


/A 


3X 


»> 1 , 


(3.13) 
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where d : is the particle diameter. Comparing it to the same term as the 
velocity gradient we find. 

| v6V 2 U i « a6K(U i -V.) 


and we can define 


K m 


2A 2v 7?U i 
3 “ ’ 3a (Uj-V,) ' 


(3.14) 


From eq. 3.14 we see that 


2 v 
3a 





« 1 


(3.15) 


is the condition to be met. Typical numbers in air are v = 0.2, a = 10?, 
(U r V i ) = 10~ 3 ; so V 2 U. « 10 3 which is not the general situation. We will 
examine the effect of a local correction considered as a constant using 
eq. 3.14 later. 

Equation 3.9 can be written as a linear first order differential 
equation using K 111 as follows: 



6 — — L + a<5K HI (U,-V.) + 6(—)% f* dT — — dx 

dt 11 71 J 

— 00 


(3.16) 


If the Basset term is dropped and K 


III 


= 1 , 


the equation becomes: 


+ a5V = a6U + 6 


dU 

dt 


(3.17) 
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The Fourier transform of the solution is 


a + iaj a, 

V = -— pU (3.18) 

at 1 6 


We follow the procedure given by Soo [2] which was developed by Tchen [4] 
and Chao [32]. Multiplying eq. 3.18 by the complex conjugates of each side 
gives 


VV 



7O 2 + ' 


(3.19) 


which is the square of the amplitude ratio for the frequency, w. For the 
particle the energy spectrum density function is introduced so that 



Hm ^ 
T^°° 2 ttT ' 


Integration with respect to w gives 



(^) 2 + 1 
a 


j < s > 2 + 1 


f (w) dco 


and 


Ie.u! 


(-) 2 

'or 


+ 1 


(“)2 

'CT 


+ 1 


(3.20) 


(3.21) 


9 


(3.22) 
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where f is the energy spectrum density function for the fluid. 

The same problem can be solved for the relative velocity and if the 
Basset term is included. Details are given by Soo. Inclusion of the 
Basset term or corrections to the drag coefficient have the same effect as 
decreasing the particle size. Therefore, the particles follow the flow 
better than would be indicated by eq. 3.21. 

Calculations have been made by Hjelmfelt and Mockros [33] for the 
amplitude at a given frequency. Results are similar to the acoustic results 
given in section II. Figure 3.1 shows the effect of particle size on the 
amplitude ratio for particles in air when all terms are included in the 
equation but only Stokes drag is considered. We have chosen cases of pipe 
flow and free jet flow where the frequency spectra is known except for the 
high frequency range for our calculations. The results given in the next 
sections are for particle-fluid relationships for the fluctuating component 
of the turbulent velocities. 

A list of the equations to be used in the subsequent discussion is given 
below so easy reference can be made to them. 

Accoustic: 



l+a 2 oj 2 
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A-II 



a = 


2r 2 p f 

9y 


Approximate Turbulent: 


T-I 


f + i 
i 7 < j > 2 + 1 


f(w) dto 


T-II 


V 2 r 

=r° I 


-(^) 2 (|) 2 


r 


o 7 <!> 2 + ' 


f (w)dto 


Turbulent with Basset force. 


T-I 1 1 


_ r°° q ^ 1 ^ 

— 2 J 

U J o 


f (oj)dw 


T-IV 


v 2 r r* ^ R (1) 

Z = / — Toy f(w)du) 

,,2 J o 


n O ) = ( W ) 2 + (“) 3/2 + 3 (— ) 

a 'a' v cr 


V5- (^) % + 1 
a 
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( 2 ) = 1 (— )2 + & (“) 3/2 + 3 (— ) + S (~) h + 1 ^ 

^2 'a' 6 a' a a 

q 0) = (L~A)2 
W R 1 5 ' V 

Turbulent with Basset force and drag coefficient (24/N^ e )K. 

fr>\ / o \ (1\ 

T-V Replace the last two terms in £r ' and ' with 


& K ) h + K 2 . 

'a' 

All velocities in the above are fluctuations from the mean. 

Figure 3.1 shows the amplitude function as a function of 

frequency and particle size for air as the fluid. Errors of over 1% exist 
when the frequency is over 6 kH 2 f° r a one m ' icron particle and for frequencies 
considerably less as the particle size increases. The effect of changing 
the parameter K is shown in Figure 3.2. As K is increased the drag increases 
and the amplitude function increases. This means that the relative velocity 
lag is reduced. 

Fortunately in a turbulent field most of the energy is concentrated 
in low frequency oscillations. That is f(co) drops off rapidly as frequency 
is increased. We examine the cases of the free jet and the center of fully 
developed pipe flow in the following subsections. The assumptions that the 
particle and fluid diffusivities are the same (long diffusion time) and no 
particle-particle interactions are retained. 

A. Application of Single Particle Analysis to Turbulent Jets 

Based on J. C. Laurence's [34] measurements in the mixing region of a 

free jet, the spectral density function can be represented by f . f 

max w 
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and a decay slope of -2 on a log-log plot of f(w) vs w as shown in Fig. 3.3. 
Using Equation T-IV extensive calculations showed the results were sensi- 
tive only to f when particle diameter and density were held constant. 

The density function decay as the inverse square cannot exist beyond 
some upper limit frequency. Most references on turbulence give a change 
from a minus 5/3 slope to a minus 7 slope for isotropic turbulence at some 
high frequency. That is at some frequency of turbulence all the energy 
should be dissipated by the viscous forces. If this frequency is taken 
as 20 kHz and we neglect everything beyond 20 kHz, calculations can be made 
giving mean square differences between particles and fluid using the equations 
listed. For a one micron diameter particle with a density of one gram per 
cubic centimeter in an air jet mixing region with f u = 2000 Hz the results 
are: 


Equation 

Result 

A- 1 

.991 

T-I 

1.000 

T-III 

.990 

A- 1 1 

.,009 

T-II 

1 009 

T-IV 

-.008 

Further calculations 

using formula T-IV are summarized in Figs 


and 3.5, Figure 3.4 shows the effect of particle density for different f 
and particle size. Figure 3.5 shows the effect of particle size when 
particle density is held constant at one gram per cubic centimeter. 


I 
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In each case the frequency spectrum density function was taken from 
Laurence's equation. 


E ( f max’ fu ’“) 


0 * fe))(fr) 2 

(4f*)[l +(^) 2 ][k, + (|y) 2 ] 


If \ 2 
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^max' 

1 2 
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iirfu ^ 



The integration was carried out using 50 points and the trapezoidal rule 
between 0 and 20 kHz. The f(u>) function in equations is the above E 
function divided by 

20 kHz 

E da) to normalize the function. 

. _2 
A simpler analysis can be made using E a linear function of a j . Figures 

3.6 and 3.7 show the ratios of particle mean or relative velocity fluctua- 
tions to the fluid fluctuation for selected particle sizes from one micron 
diameter to 200 micron diameter. Air is taken as the fluid at room temper- 
ature. 

The effects of non Stokes drag, rarefaction, compressibility and 
inertia can be included in the equation of motion by letting 
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(1 + 0.15N R °' 687 )(1 + exp (-0. 427/M 4 ' 63 - 3.0/N R 8 ' 88 )) 

1 + /- (A + B exp (-1.25N Re /M)) 

Rg 

< • 

where the Reynolds number, N Rg , and the Mach Number, M, are for the 
difference velocity between the fluid and particle. The constants A and 
B are 3.82 and 1.28 respectively at Mach numbers sufficiently high, but 
may be lower for lower Mach numbers. The equation above is derived by 
Carlson and Hoglund [27]. 

For a one micron diameter particle with a density of one gram per 
cubic centimeter moving in air at room temperature, the following results 
were obtained: 


Mach No. 

Mean square 

Mean square 

Particle 

Particle 

of fluid 

ratio of 

relative velo- 

Reynolds 

Mach 


velocities 

city to fluid 
velocity 

Number 

Number 

0.1 

.9878 

.0096* 

.218 

.0098 

0.3 

.9889 

.0087* 

.623 

.0281 

0.5 

.9897 

.0082* 

1.003 

.0452 

1.0 

.9909 

.0072 

1.883 

.0848 

2.0 

.9924 

.0061 

3.454 

.1556 


In the table the complete equation was used although the rarefaction 
correction in the denominator is not valid for the starred points. In 
general for this particle size at a cut-off frequency of two KHz as used 
in the calculations there is no significant change. All of the conditions, 
however, fall in the slip flow regime and not the continuum regime. 

B.~ Single Particle in Pipe Flow 

A similar analysis can be made for fully developed turbulent pipe 
flow. For the frequency spectra Laufer’s data [35] can be used. The upper 
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cut off frequency is now a function of Reynolds number or average velocity. 
If the minus second power decay of the frequency spectra can be used, the 
curves for the free jet would also be valid for the drag corrections. From 
Laufer's curve at the centerline of the one dimensional spectra 


0.28 U • = fu 
X-) max 

where fu is in radians per second and U in cm. /sec. Laufers curve 

X*| maX 

is for a Reynolds number of 5 x 10^, however, and may not be representative 
of lower Reynolds numbers. 


C. Experimental Studies of Turbulent Velocities of Small Particles 

A summary of the few experimental measurements on turbulence inten- 
sities in particle-fluid flows is given by Reddy and Pei [36]. Other 
studies discussed in Part IV of this report are indirect in that measure- 
ments were made on concentration fluctuations and not velocities. The 
velocity studies were made photographically using particles greater than 
100 microns. Reddy and Pei conclude that "the particle turbulence compo- 
nents can be correlated principally with the fluid turbulence." It is 
difficult to justify this statement based on the calculations in this 
report. Only the slowly varying large scale motions follow the flow and 
therefore the turbulence intensities of the particle and the fluid check 
closely. However, it is impossible to measure the high frequency turbulence 
components by measuring velocities of particles larger than one micron. 

A study in liquids using 0.5 micron particles was recently reported 
by Corino and Brodkey [37]. They were successful in measuring particle 
velocity fluctuations near the wall of a pipe at Reynolds numbers up to 
50,000. Velocity gradients as high as 10^ sec.“^ could be determined. 
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Soo [2] reviews other experiments but the interaction of the solid 
particles with the wall has been a problem. Min [38] indicates that the 
spread of velocities after collision with the wall must be considered and 
not confused with the turbulence of the main stream. In section I of this 
report the time decay was seen to influence the relative mean velocity. 
Measurements of turbulence from velocities of large particles would con- 
sist of large effects due to these adjustments following collisions. 



IV. OTHER CONSIDERATIONS IN PARTICLE-FLUID FLOWS 


Two assumptions of the previous section are not substantiated by 
experiments. These are the assumption that particle and fluid diffusivities 
are the same, and the neglect of particle-particle collisions. Extensive 
theoretical and experimental investigations have appeared in the literature 
on both problems, but it is not possible to predict or measure the effects 
in most real situations. 

Soo [2] derives the eddy diffusivity for the particle from the mean 
square displacement. If the mean square displacement is related to the 
Lagrangian velocity auto correction coefficient, R(t), the result is 

D p (t) = V 2 P R (T)dx (4.1) 

o 


The fluid eddy diffusivity similarly is 



R(t) dx 


(4.2) 


Friedlander [29] explains that R(x) is Lagrangian only for small particles 
and may be between the Eulerian and Lagrangian correlations for the real 
case. The ratios of particle to fluid eddy diffusivity become 


D "2 


for short times 



for long times 


and 
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Both Soo and Friedlander give results for an assumed Lagrangian coefficient. 
In the previous section the assumption was made that Dp/D = 1, but studies 
by Lumley [28] and others show that this is not true. To account for the 
difference we must relax the assumption that a particle remains always with 
the same fluid element. Soo gives a review of calculations to correct the 
eddy diffusivity for some practical cases. If the particles are below one 
micron in diameter, however, the assumptions of section III appear to be 
valid when sufficient time is given for the particle-fluid system to 
equilibrate. Collisions with walls or other particles, velocity gradients 
or coalescence will affect this time period. 

Some of the first work in the field of particle-particle interactions 
was published by Smoluchowski in 1916. He developed an infinite set of 
simultaneous differential equations to describe the interactions and growth 
of floes. His results may be written in condensed form as 


dN k 

dt 


rk-1 


2 dU 

3 dZ 


L ww 1 


Li=l 

j=k-i 


2N k I MW 


i=l 


(4.3) 


where dN/dt is the rate of change of the number of particles of the indicated 
size, dU/dZ is the velocity gradient, N is the number of particles of an 
indicated size, and R is the radius of the particles. It is not possible 
to solve this set of equations in a closed form for real cases; however, 
with the advent of large computer systems approximate solutions have been 
made. [ 39 , 4(3* 
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A program is available from R. S. Gremmell [41] in which the 
growth, breakup, and steady state size of floes in laminar flow are, 
simulated as a function of the maximum floe size, the breakup floe sizes, 
the velocity gradient, and the initial concentration of the suspension. 

This type of simulation is very helpful if the necessary variables are 
known for a particular system. Such simulations are, of course, only 
useful if the particles tend to floe or coagulate. 

The availability of sophisticated equipment has made experimental 
studies of single particle interactions possible. In the study of rigid 
spheres, Mason and others [42] discovered a repeatable pattern of approach, 
interaction, and separation. The spheres approach each other rectil inearly 
in a path parallel to the x-axis (i.e. axis parallel to the flow). Although 
they never come into true physical contact the two approaching spheres 
form a doublet which rotates with a constant angular velocity of 



where G is the velocity gradient in the fluid. The doublet rotates about 
an axis perpendicular to the planes of shear, and in the case of rigid 
spheres the doublet only rotates 180°. The separation is a mirror image 
of the approach because of the 180° rotation before separation. Deform- 
able spheres do not separate after the first half cycle but may rotate 
for several full cycles before coalescing or separating. They will 
coalesce only if the van der Waals forces are strong enough to displace 
the intervening liquid. [ 43 ] 

The tubular pinch effect is a fluid induced particle disturbance 
found only in tubular Poiseuille flow and not in Couette flow. The 
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tubular pinch effect is essentially the tendency for rigid or deformable 
bodies to migrate under certain conditions toward a radial position of 
distance r from the center of the tube. The particles closer to the 
center of the tube migrate inward. This effect is apparently not en- 
countered in turbulent flows although the forces causing the phenomenum 
may be important in the analysis given in the previous section. 

Karnis and Mason [44] have obtained an expression for the radial 
position, r, of a liquid drop as a function of time, t, and indicated 
that a rigid sphere follows this behavior very closely. Their result 
is given as 


R 


-{r 

0 0 


k 2 b 4 n o t 


yR, 


0 I 


33(54p^+102p+54) (19p+16) + A (44) 
4480(p+l) 3 


where R q is the radius of the tube, k is given as 4Q/(R 4 ), Q being the 
volumetric flow rate, b is the original radius of the undeformed drop, n 0 
is the viscosity of the suspending fluid, y is the interfacial tension, p 
is the viscosity ratio of the suspended to the suspending phase, and 




) 


where r Q is the initial position of the drop. The actual migration rate 
in some cases has been found to be lower than that predicted by equation 

4.4. 

Takano, Goldsmith, and Mason [45] investigated the effect of a sine 
wave fluid disturbance superimposed on the steady state flow. Although 
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their work is not complete, they have come to some qualitative conclusions. 
The migration rates and equilibrium positions for a given suspension 

t 

depend on a dimensionless parameter, a, as well as the flow rate. The 
parameter 


where R q and n Q are the same as before, w is the angular frequency of 
oscillation, and p Q is the density of the suspending fluid. This parameter 
was used in the analysis of both rigid and deformable spheres. 

In experiments made with polystyrene spheres, the initial stages 
of oscillatory flow, a<5, revealed an increased migration rate over that 
of steady state Poiseuille flow. The particles collided in approximately 
the same manner as observed in the former particle interaction studies in 
steady state laminar flow. As a became larger than ten, more than one 
migration or equilibrium position was noted. Some of these positions, 
defined as radii where dr/dt for the particle equaled zero, were found to 
be quasi-stable. The number of these positions increased with increasing 
a, and it would be expected that at sufficiently high a the number of these 
quasi-stable positions would destroy any noticeable migration caused by 
the tubular pinch effect. 

In a study concerned with the agglomeration of aerosols induced by 
high energy sound waves St. Clair [46] listed three fluid effects on 
particle movements. These were convibrational , hydrodynamic and radiation 


pressure. 



Another analysis by Shaffman and Turner [47] concerned the collision 
of liquid droplets in turbulent clouds. Assuming isotropic turbulence, the 
collisions of the drops are attributed to two main sources: (1) collisions 
due to the motions of the droplets in (within) the air, (2) collisions due 
to motions of the droplets relative to the air (i.e. inertial effects). 

The two effects are discussed separately; the two effects may be calculated 
separately and then added together. 

The collisions caused by the drops moving with the fluid are discussed 
assuming spherical droplets with a sphere of influence, R, equal to the sum 
of the radii of two particles, r-j + r 2 - With a radial particle velocity 
of w(r), some fluid relations can be used to transform w(r) into fluid 
properties. After the transformation, the number of collisions, N, bet- 
ween two population concentrations, n-j and n^, is 

N = 1.3 n^ 2 (r 1 + r 2 ) (£) , (4.5) 

where e is the rate of energy dissipation per unit mass, and y is the 
kinematic viscosity. 

In order to obtain the total effect of the collisions, the collision 
rate is substituted into the Smoluchowski equation. 

In dealing with the second part of the problem, that of droplet 
motions relative to the air as causes of collisions, a probability function 
P(w)for the relative velocity of the particles is assumed and the collision 


rate thus obtained is 
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By introducing a collision efficiency a quantitative prediction of droplet 
formation in turbulent clouds can be obtained. 

< 

Another theoretical development of turbulent coagulation was made 
by H. C. Frish [48]. Homogenous and isotropic turbulence was considered 
in a fluid macroscopically at rest. In order to describe such turbulence, 
a factor of turbulent diffusion, n* is defined as 

N*(t) = v 2 / R h (a)da + D (4.7) 

o 

where D is the molecular diffusion coefficient, R h is the Lagrangian 
correlation coefficient of the field of turbulence and v is the variance 
of the turbulent velocities. This definition is augmented with the assump 
tion that the probability density for the particle is Gaussian and which 
gives 

N T, 2 (t) = N*(t) + N*(t) (4.8) 

He next makes the common assumption that a particle fixed in space 
has a sphere of influence of R-j 2 = r i + r 2 > where r-j and r ^ are the res- 
pective radii of two colliding particles. He considers the collision 
frequency to be one. Using the above conditions an expression for the 
rate of formation of a particle of the size i+k can be written. Then 
the concentration distribution can be calculated. 

Becker, Hottel, and Williams, optically, measured the concentration 
fluctuations in a particle-fluid nozzle flow [49]. Although they assumed 
the velocity and concentration fluctuations were alike, the results showed 
the characteristic difference in the diffusion broadening. It is apparent 
that little data are available to properly check the differences between 
particle and fluid velocity fluctuations at high rates of change. 



V. CONCLUSIONS 


Small particles can be used to measure the fluid velocity motions 
except under conditions of high frequency fluctuations.. At high velocities 
in supersonic flow and for detailed analyses of turbulence, all particle- 
fluid-instrument interactions must be considered. Both more theoretical 
correlations for these conditions and experimental measurements under more 
ideal conditions are needed to determine the exact limits of the Laser- 
Doppler Flowmeter and the particle lag corrections. 
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